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The two-photon width of \c2 ( P2) state of charmonium has been measured using 14.4 ftT 1 of 
e + e~ data taken at y^i = 9.46 — 11.30 GeV with the CLEO III detector. The 77-fusion reaction 
studied is e + e~ — > e + e~77, 77 — » Xc2 — > ~yJ/ip — » je + e~ (fi + n~). We measure Y 11 {xc2)B{xc2 — > 
jJ/tjj)B(J/ijj -> e+e" + fi + H~) = 13-2 ± 1.4(stat) ± l.l(syst) eV, and obtain T 77 (x C 2) = 559 ± 
57(stat) ± 48(syst) ± 36(br) eV. This result is in excellent agreement with the result of 77-fusion 
measurement by Belle and is consistent with that of the pp — > Xc2 — > 77 measurement, when they 
are both reevaluated using the recent CLEO result for the radiative decay Xc2 — > yJ/ip- 

PACS numbers: 13.20.Gd, 13.40.Hq, 14.40.Gx 

The P-wave states of charmonium ( 3 Pj, 1 Pi) have provided valuable information about the qq interaction 
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and QCD. The two-photon decays of the positive C- 
parity states ( 3 Pj) are particularly interesting because 
at lowest order the two-photon decay of charmonium is 
a pure QED process akin to the two-photon decay of 
positronium. Their study can shed light on higher order 
relativistic and QCD radiative corrections. 

The measurement of the two-photon width of Xc2, the 
3 P 2 state of charmonium, r 77 (x C 2) = T(Xc2 — > 77), has 
a very checkered history, with large differences in results 
from measurements using different techniques. The pre- 
1992 77-fusion measurements of r 77 (x C 2) were inconclu- 
sive, with most of them only establishing upper limits 
of several keV. In 1993, the E760 experiment at Fermi- 
lab reported the result from their pp —* Xs2 ~ * 77 mea- 
surement, r 77 (xc2) = 320 ± 80 ± 50 eV a factor of 
more than 3 smaller than the smallest limit established 
by the 77-fusion measurements. The 77-fusion experi- 
ments continue to report much larger valuers of r 77 (x C 2) 
than the pp experiments, with the result of the 77-fusion 
measurement from Belle 2], r 77 (x c2 ) = 850 ± 127 eV, 
which is still three times larger than the latest pp mea- 
surement of Fermilab E835 0, r 77 (x c2 ) = 270±59 eV. It 
is this continuing discrepancy between the present good- 
statistics measurements that has motivated the investi- 
gation reported here. 

In this investigation we report on a measurement of 
r 77 (x C 2) by the study of the 77-fusion reaction 

e + e _ -> e+e~(77) , 77 -> \c2 -> jJ/ip -» (1) 

The data sample used for the analysis was taken at 
the Cornell Electron Storage Ring (CESR) with the de- 
tector in the CLEO III configuration The detector 
provides 93% coverage of solid angle for charged and neu- 
tral particle identification. The detector components im- 
portant for this analysis are the drift chamber (DR) and 
CsI(Tl) crystal calorimeter (CC). The DR and CC are op- 
erated within a 1.5 T magnetic field produced by a super- 
conducting solenoid. The DR detects charged particles 
and measures their momenta. The CC allows measure- 
ments of electromagnetic showers with energy resolution 
a(E)/E = 2.3 - 2.7% for E 1 = 0.3 - 0.6 GeV. 

The data consist of a 14.4 fb^ 1 sample of e + e~ colli- 
sions at or near the energies of T(15 — 55) resonances 
and around the AjAft threshold in the range of center-of- 
mass energies yfs = 9.46 — 11.30 GeV. The data sample 
sizes are given in Table I. 

The two-photon partial width r 77 (x C 2) was measured 
in untagged 77-fusion reaction of Eq. (1). Events with 
r ye + e~ or r yp + fj,~ in the final state were selected. The se- 
lected events are required to have two charged tracks and 
zero net charge. All charged particles were required to 
lie within the drift chamber volume and satisfy standard 
requirements for track quality and distance of closest ap- 
proach to the interaction point. 

The photon produced in the decay Xc2 — ► ~fl + l~ typi- 
cally has an energy E 1 k, 0.46 GeV. The selected events 



TABLE I: Data used in the present analysis. Average values 
of y/s and corresponding luminosities are listed. 



Data 


T(15) T(25) T(35) T(45) T(55) AjAf, 


C (ft)" 1 ) 
(GeV) 


1.399 1.766 1.598 8.566 0.416 0.688 
9.458 10.018 10.356 10.566 10.868 11.296 



are required to have only one electromagnetic shower 
with an energy 0.3 < E 1 < 0.6 GeV, and to be isolated 
from the nearest charged track by an angle > 20°. The 
total energy of remaining electromagnetic showers in the 
event £ ; to t(neut) was required to be < 0.3 GeV. 

The total energy of the system, _E to t(7^ + ' - ), defined 
as the energy sum of the lepton pair and the candidate 
photon, was required to be less than 5 GeV. This cut has 
an efficiency of ~ 96% and removes all background that 
arises when "0(25) is produced via initial state radiation, 
decays to jjJ/tp, and one photon is not detected. 

Untagged 77-fusion events are characterized with 
small transverse momentum; therefore p^- ot ( r yl + l~) < 
0.15 GeV/c was required. 

Lepton pairs of low transverse momentum may also be 
directly produced by 77-photon fusion. These constitute 
a background that is removed by rejecting lepton pairs 
withp^Z+r) < 0.1 GeV/c. 

To identify two charged tracks as electrons or muons, 
the E/p variable was used, where E is the energy de- 
termined from the calorimeter and p is the momentum 
determined from track reconstruction. For both muons 
< E/p < 0.3 is required, and for both electrons 

0. 85 < E/p < 1.15 is required. If a photon of energy 
larger than 0.03 GeV is present within a 5° angle cone 
around the lepton direction, it is assumed to be the re- 
sult of bremsstrahlung and its momentum is added to the 
momentum of the track. 

The signal Monte Carlo (MC) sample for untagged 77 
fusion production of Xc2 resonance was generated using 
the 77 fusion formalism from Budnev et al. MC sam- 
ples were produced for each value of yfs listed in Table 

1. For the calculation of the overall event selection ef- 
ficiencies, MC samples were weighted according to the 
luminosity of each data set. 

According to Budnev et al. 0, when the photons are 
transversely polarized, the untagged 77 production cross 
section is related to the two-photon cross section by 

dcr e + e -^ e + e - Xc2 = c?£ 77 r \W 2 )a^^ Xc2 , 

where C'EI' is the 77 luminosity function and W is the 
two-photon invariant mass. The implementation of the 
above formalism has been discussed in detail by Dominick 
et al. ;6|, and we follow it. We calculate cr(x c2 )/r 77 (x C 2) 
for each value of yfs assuming that Xc2 production in the 
fusion of two transverse photons is significant only be in 
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the helicity 2 state 0], and that the radiative transition 
Xc2 — > jJ/ip is pure El. We assume that the intermediate 
vector meson in the Budnev formalism is a J/ip, and we 
implement the proper angular distribution ||J in calculat- 
ing efficiencies. The luminosity- weighted average value of 
c(Xc2)/r 77 (x C 2) is determined to be 4.93 pb/keV. 

Good agreement is observed between the data and MC 

distributions for £? t ot(7^ + ' _ ), Ptot (7^ )' ^71 ^~ mo ~ 
menta, -Etot(neut), E/p for leptons, and the photon and 

lepton angular distributions. The E/p and p^ t (7Z + Z~) 

distributions are shown in Fig. 1, and the lepton and 

photon angular distributions are shown in Fig. 2. 

For the data at each a/s, efficiencies of all event selec- 
tion requirements were determined from signal MC sim- 
ulations, and the averages are listed in Table II. 

A two dimensional plot of the lepton pair masses 
(e + e~ + versus the mass difference AM — 

M(jl + l~) ~ M(l+l~) is shown in Fig. 3. It shows a 
clear enhancement at the mass of J/ip. A cut M(l + l~) = 
M (J/ip) ±30 MeV was therefore used. The resulting dis- 
tributions of AM are shown in Fig. 4 for (a) e + e~, (b) 
and (c) e + e~ plus 

Three different methods, all using the back- 
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FIG. 1: Distributions of E/p (left) and p tot (right) in data 
(points) and in the signal MC (histograms). Vertical dashed 
lines indicate the E/p cuts for electrons and muons. 
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FIG. 2: Comparison of the background subtracted data 
(points) and the signal MC (histograms) distributions of 
jcos0 7 7 j, jcosO; + j and A(p. 7 7 is the polar angle of the 
photon in the rest frame, O l l+ is the polar angle of the 

positive lepton in the rest frame, and Atp is the azimuthal 
angle difference between the momenta of the two leptons in 
the laboratory frame. 



TABLE II: Efficiencies of the different event selection criteria. 



Selection cut 


e + e 






Channel (%) Channel (%) 


NY rViaro-p^l — 9 

1 > l i^lldl tit; j — _■ 


68.9 


70.8 


Total Charge=0 


98.7 


98.7 


Only one *y with 






0.3 < E 1 < 0.6 GeV 


52.8 


53.7 


T .PTiTr^n / T) 


92.4 


98.3 


Etotijl+r) < 5 GeV 


96.1 


95.3 


Etot(neut) < 0.3 GeV 


99.0 


99.1 


pL(l + l~) >0.1 GeV/c 


99.0 


98.9 


pL(fl + l~) < 0.15 GeV/c 


62.1 


62.4 


M(l+r) = M(J/ip) ± 30 MeV 


81.9 


93.0 


Trigger 


97.5 


85.7 


Overall efficiencies 


15.5 


17.1 



ground shape determined from the J/ip sideband re- 
gion [Af(Z+Z-)=2.7-3.5 GeV, omitting M(l + l~ )=3.0-3.2 
GeV], were used. Fits using the Crystal Ball line shape 
(which consists of a Gaussian with a low energy tail) |3] , 
signal MC peak shape, and simple counts in the region 
AM=0.42-0.49 GeV led to yields and efficiencies that 
differ by less than ±2%. 

The observed yield, N a b s , leads to 



r 77 (x c2 ) x B(xc2 -» ji + r) 



N 



bs 



eC{cF{xc2)/T 11 (xc2)Y 



where e is the total efficiency, C is the total luminos- 
ity of the data used, c(Xc2)/r 77 (xc2) = 4.93 pb/keV as 
determined earlier, and B(x C 2 —* ll + l~) = Bi(Xc2 — * 
jJ/iP)B 2 (J/ip -» 1+1-). 



3. 
+ 



3.6 
3.5 
3.4 
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FIG. 3: Scatter plot of the AM = M(jl+l~) - M(l + l~) 
with respect to the two lepton effective mass in data. Top 
and bottom horizontal shaded bands are the areas defined as 
the J/ip sideband regions, and the middle band is the area 
defined as the J/ip signal region. 
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TABLE III: Average of results for the three signal count 
extraction methods. 



0.4 0.5 0.6 

AM (GeV) 

FIG. 4: AM = M(jl + r) - M(l + r) mass difference distri- 
butions for e + e~ channel (a), /x + \i~ channel (b), and the sum 
(c). The solid line curves are results of fit to the data points 
using the background indicated by the dashed line curves. 



In Table III we present our results, which are averages 
of the results for the three different signal yield extraction 
methods. We present the results for e + e~ and \x + \x~ 
separately, and for their sum. Our directly determined 
result for the sum is 



i + r N obs r 77 ( Xc 2)z?(xc2 ~*ji + r 


") r 77 (xc2) 


(eV) 


(eV) 


e+e" 68±11 6.4±1.0 


544±87 


79±11 6.8±0.9 


571±76 


Total 147±15 13.2±1.4 


559±57 


TABLE IV: Sources of systematic uncertainties. 


Source Systematic 


uncertainty (%) 


integrated luminosity, C 


±3.0 


trigger efficiency 


±3.0 


signal yield extraction 


±1.3 


J/ip line shape modeling 


±1.6 


photon resolution modeling 


±1.3 


event selection 


±4.8 


tracking 


±2.0 


photon finding 


±2.0 


J/ip (versus p, <f)) in 77 


±3.0 


pure El (versus El + 10% M2) 


±3.0 


overall 


±8.6 



anticipated by Belle @ Both Belle @ and E835 used 
the 2000 PDG value of B{ Xc 2 -> 7J/VO = (13-5 ± 1.1)%. 
As shown in Table V, when these results are reevaluated 
using B(xc2 — * jJ/tp) — (19.9 ± 1.3)%, as recently mea- 
sured by CLEO [9j, the Belle result [2| comes into com- 
lete agreement with ours, and even the latest pp result 
becomes statistically consistent with ours. 
Many theoretical predictions based on potential model 
calculations exist in the literature. In an early relativis- 
tic calculation, Barnes predicted r 77 (x c2 ) = 560 eV. 
In calculations including both relativistic and one-loop 
QCD radiative corrections, Gupta et al. 0] and Ebert 



r 77 (Xc2)#(x C 2 7(e + e~ ± /i + /0) = 13.2 ± 1.4 eV. et al predict r 77 (x c2 ) = 570 eV and 500 eV, respec 



We use B\{xc2 —* jJ/ip) — (19.9 ± 1.3)% as measured by 
CLEO 0, and [Bn{J/ip -> e+e~) = B 2 (J/ip -> = 
5.95 ± 0.06)% [Tj, to obtain r 77 (x c2 ) listed in Table III. 
Various sources of systematic uncertainty were studied. 
These are listed in Table IV, and are combined in quadra- 
ture to give a total systematic error of ±8.6%. Thus, our 
final result is 

r 77 (x c2 ) = 559 ± 57(stat) ± 48(syst) ± 36(br) eV. 

Table V shows a compilation of the published re- 
sults of Belle @, E835 and CLEO 13], together 
with our result. We find that a large part of the dis- 
crepancy between the earlier 77-fusion results and the 
PP ~ * Xc2 —> 77 results arises from the use of the old 
values of B(xc2 ^ -fJ/ip). This possibility was indeed 



tively. All three predictions are in agreement with our 
result. 

An estimate of the strong coupling constant, as{m c ) 
can be obtained by comparing r 77 (x C 2) with r 9S (x C 2)- 
With the known first order QCD radiative corrections 
for the two widths [13, the pQCD prediction is 



r 77 (x c2 ) 

T gg(Xc2) 



80^ 
9a? 



1 



(2) 



The hadronic width, r gs (x c2 ) = r to t(Xc 2 ) X B(xc2 -> 
gg) = T tot ( X c2)x[l-B(xc2 ^ 1J/1P)] = 1.55±0.11MeV, 
obtained by using r to t(Xc 2 ) = 1.94 ± 0.13 MeV [ill and 
B(xc2 -> jJ/ip) = 0.199 ± 0.013 0|. Using our result for 
r 77 (Xc 2 ), we obtain r 77 (x c2 )/r s9 (x c2 ) = (361 ± 59) x 
10~ 6 . Equating this to the pQCD expression (Eq. 2) but 
not including the QCD radiative corrections in the large 
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TABLE V: Comparison of our result for F 77 (xc2) with the 
results of the two recent 77-fusion measurements and the 
Fermilab E835 pp experiment. The second column gives 
the results as published and the third column gives the re- 
sult after reevaluation using the CLEO measured values for 
Bi(Xc2 -> 7J/VO and B 2 (J/ip -> l + l~) 0,E|. Also, the aver- 
age Fermilab measured value of r to t(Xc2) HU is used to recal- 
culate the E835 result 0, and PDG2004 value of B{xc2 -> 4tt) 
is used to recalculate the CLEO result [l3|. 



Experiment [Ref.] 
Quantity Measured 



F 77 (Xc 2 ) (eV) r 77 ( X c2) (eV) 
(as published) (as reevaluated) 



Present: 77 — > \ c :i 

jl+l-) 



F 77 (Xc2)S(Xc2 



559(57)(48)(36) 



Belle [2]: 77 - 

F 77 (Xc2)S(Xc2- 



-» Xc2 

>7' + r) 



850(80)(70)(70) 570(55) (46) (37) 



CLEO [13]: 77 

T 77 (Xc2)S(Xc2- 



* Xc2 
>4tt) 



530(150) (60) (220) 432(122) (54) (61) 



E835 [3]: PP~*Xc2 
(Xc2->77)/(Xc2->7-W0 



270(49)(33) 



384(69) (47) 



parentheses, gives a s (m c ) = 0.36 ± 0.03. The pQCD 
expression with the QCD radiative corrections (Eq. 2) 
leads to the value a s (m c ) = 0.29 ± 0.02. 

We gratefully acknowledge the effort of the CESR staff 
in providing us with excellent luminosity and running 
conditions. This work was supported by the National 
Science Foundation and the U.S. Department of Energy. 
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